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Energy and exergy assessments are reported of a novel trigeneration system based on a solid oxide fuel 
cell (SOFC), for steady-state operation and using a zero-dimensional approach. The trigeneration system 
also includes a generator-absorber heat exchanger for cooling and a heat exchanger for the heating 
process. The influences of two significant SOFC parameters (current density and inlet flow temperature) 
on several variables are investigated. The results show that the energy efficiency is a minimum of 33% 
higher when using the trigeneration system compared with the SOFC power cycle. In addition, the 
maximum energy efficiencies are found to be 79% for the trigeneration system, 69% for the heating 
cogeneration, 58% for cooling cogeneration and 46% for electricity production. Moreover, the highest 
trigeneration exergy efficiency is almost 47% under the given conditions. It is also shown that, as SOFC 
current density increases, the exergy efficiencies decrease for the power cycle, cooling cogeneration, 
heating cogeneration and trigeneration. As current density increases, the trigeneration energy and exergy 
efficiencies decrease, and an optimal current density is observed to exist at which the net electrical power 
is a maximum. As SOFC inlet flow temperature increases, the trigeneration energy and exergy efficiencies 
and net electrical power increase to a peak and then decrease. The main exergy destructions occur in the 
air heat exchanger, the SOFC and the afterburner. 

© 2014 Published by Elsevier Ltd. 


1. Introduction 

Accessibility of energy resources and global warming are two 
important long-term problems. Global energy demand has steadily 
risen over time, notwithstanding the limited availability of non¬ 
renewable fossil fuels. As a consequence, more efficient and clean 
energy systems are being sought. The energy efficiency of conven¬ 
tional power plants based on single prime mover is usually less 
than 40%, with most of the energy loss associated with waste heat. 
Combining cooling and heating subsystems in a conventional plant 
can increase the system efficiency considerably, and energy effi¬ 
ciencies can reach 80% through trigeneration [1,2], the simulta¬ 
neous production of cooling, heating and electricity, often from 
one energy source. Trigeneration, also referred to as a combined 
cooling, heating and power (CCHP), offers a number of potential 
benefits: 
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• Increased efficiency. 

• Reduced fuel and energy costs. 

• Significantly reduced greenhouse gas emissions. 

• Shortened transmission lines. 

• Fewer distribution units. 

• Increased power supply reliability. 

The primary driver is an important part of a trigeneration 
system, making its selection important. The main primary drivers 
are gas turbines, steam turbines, external combustion engines, 
internal combustion engines, micro-turbines, and fuel cells. Among 
these, trigeneration based on fuel cells often achieves higher 
energy efficiency, in part because the fuel cell avoids Carnot 
efficiency limitations. 

The solid oxide fuel cell (SOFC) is considered an important 
emerging technology. It can produce electricity directly from a fuel, 
such as methane, and an oxidant, such as air. It operates at high 
temperatures and generates much waste heat, so it is commonly 
coupled with a gas turbine (GT) or an organic Rankine cycle 
(ORC) as a bottom cycle to improve the overall efficiency by 
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Nomenclature 



a 

activity factor of species in electrochemical reaction 

Greek letters 

A, 

active surface area, m 2 

>1 

energy efficiency 

^aeff 

effective gaseous diffusivity through anode, m 2 /s 


exergy efficiency 

^ceff 

-ex? 1 ’ 0 

effective gaseous diffusivity through cathode, m 2 /s 
standard chemical exergy of species, kj/mol 

p 

electrical resistivity of cell components 

Ex 

exergy rate, kW 

Subscripts 

£x ch 

rate of chemical exergy, kW 

0 

atmospheric condition 

£x ph 

rate of physical exergy, kW 

a 

anode 

f 

Faraday constant, C/mol 

abs 

absorber 

Ag° 

change in molar Gibbs free energy, J/mol 

ac 

AC current, air compressor 

GAX 

generator absorber heat exchanger 

act 

activation 

h 

molar enthalpy, J/mol 

AHE 

air heat exchanger 

I 

current, A 

c 

cathode 

i 

exergy destruction rate, kW 

C0g,c 

cooling cogeneration 

i 

current density, A/m 2 

cog,h 

heating cogeneration 

Jas 

anode-limiting current density, A/m 2 

cone 

concentration 

jcs 

cathode-limiting current density, A/m 2 

cond 

condenser 

J oa 

exchange current density of anode, A/m 2 

cv 

control volume 

Joe 

exchange current density of cathode, A/m 2 

des 

desorber 

I< 

equilibrium constant of waste gas shift reaction 

e 

electrolyte, exit 

L 

thickness of SOFC layer, m 

el 

electrical power 

LHV 

lower heating value, kj/mol 

evap 

evaporator 

n 

molar flow rate, mol/s 

fc 

fuel compressor 

Nfc 

total number of fuel cells 

FC 

fuel cell 

Q 

heat rate, kW 

FHE 

fuel heat exchanger 

R 

universal gas constant, J/mol K 

h 

heating 

^contact 

contact resistivity, O m 2 

HR 

heat recovery 

f*el,h 

electrical to heating ratio 

in 

inlet 

^el.c 

electrical to cooling ratio 

int 

interconnect 

S 

molar entropy, J/mol K 

inv 

DC to AC inverter 

T 

temperature, K 

ohm 

ohmic 

u f 

fuel utilization ratio 

N 

Nernst 

u 0 

air utilization ratio 

tri 

trigeneration 

V 

voltage, V 

WHE 

water heat exchanger 

w 

power, kW 

wp 

water pump 

Wpc, stack 

power output of fuel cell 

_ 

molar 

*r 

extent of steam reforming reaction for methane, mol/s 


rate 

yi 

molar concentration 

0 

standard pressure 

y r 

extent of water gas shift reaction, mol/s 

ch 

chemical 

Zr 

extent of electrochemical reaction, mol/s 

ph 

physical 


recovering SOFC waste heat. Many researchers have studied the 
performance of SOFC-GT hybrid power plants [3-16], integrated 
SOFC-ORC power plants [17,18] and combined SOFC-Stirling 
hybrid plants [19J. Few analyses, however, have been reported of 
trigeneration power plants based on an SOFC. Some of the more 
significant of these are now described. 

Weber et al. [20] conducted detailed C0 2 emission and cost 
analyses of a trigeneration system based on a SOFC primary mover 
in an office building, and demonstrated that the system reduced 
C0 2 emissions by 30% at a cost increase of 70% compared with a 
conventional system. 

A trigeneration plant driven by a SOFC and a gas turbine was 
studied by Burner et al. [21]. The plant included half-, single- and 
double-effect chillers, a heat pump, an additional gas boiler, a heat 
recovery boiler, a compression chiller and cooling system. They 
considered the potential of combining a heat pump into the 
trigeneration system. The study considered energy and exergy effi¬ 
ciencies and concentrated on costs and C0 2 emissions, and demon¬ 
strated that the SOFC-GT system is attractive economically and 
environmentally. The study utilized multi-objective optimization 


based on a multi-criteria evolutionary algorithm, and incorporated 
the following optimization objectives: annual total cost of trigener¬ 
ation and annual C0 2 emission rates. 

Liu et al. [22] proposed a trigeneration plant by combining an 
internal-reforming solid oxide fuel cell (1RSOFC) with a zeolite/ 
water adsorption chiller, and analyzed numerically the perfor¬ 
mance of the system under various operating conditions and 
parameters. 

Yu et al. [23] analyzed a total energy system incorporating a 
SOFC and an absorption chiller driven by exhaust gas to provide 
electricity, cooling and/or heating simultaneously. They found that, 
as fuel utilization coefficient varies, both the electrical efficiency 
and total efficiency of the system achieved maximum values. 

Al-Sulaiman et al. [24] carried out an energy analysis of a tri¬ 
generation system based on a SOFC and an ORC. The system also 
consisted of a heating process and a single-effect absorption chil¬ 
ler. The analysis demonstrated that, compared to only a power 
plant, trigeneration increases the efficiency by at least 22%. They 
also found that the exergy efficiency increased by 3-25% when 
trigeneration was applied [25]. 
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Wang et al. [26] conducted energy and exergy analyses of a 
methane-driven trigeneration system based on a SOFC. The system 
included an IRSOFC, an afterburner, a gas turbine, heat exchangers, 
compressors and a Kalina cycle. A parametric study indicated that, 
as air flow rate increases, the SOFC electrical efficiency, overall 
electrical efficiency and trigeneration exergy efficiency all increase. 

A trigeneration system integrating a SOFC with a gas turbine 
was investigated with energy analysis by Shaolin Ma et al. [27], 
for various operating conditions and parameters. 

In the present study, energy and exergy analyses are carried out 
of a proposed novel trigeneration plant based on a SOFC, under 
steady-state operation and using a zero-dimensional approach. In 
the zero dimensional model, it is assumed that the fluid state is 
homogenous in the control volumes. This trigeneration system 
consists of an SOFC, a generator-absorber heat exchanger (GAX) 
for cooling and a heat exchanger for heating. The influence is inves¬ 
tigated on six variables (cell voltage, voltage loss, electrical-power, 
heating cogeneration, cooling cogeneration and trigeneration) of 
variations in such parameters as SOFC current density, SOFC inlet 
flow temperature and fuel utilization factor. 

2. System description 

A schematic diagram is shown in Fig. 1 of the proposed novel 
trigeneration system based on a SOFC as the primary driver. The 
system includes a SOFC, a GAX and a heating-process heat 
exchanger. The exhaust gases from SOFC system supply heat to 
the generator of the GAX system for cooling purposes. Then, resid¬ 
ual heat is used in the heat exchanger to produce hot water. 


3. Assumptions 

Several assumptions are employed to simplify the present inves¬ 
tigation. The analysis is carried out under thermodynamic equilib¬ 
rium and steady state conditions. In addition, changes in kinetic and 
potential energy are neglected since they are relatively minor. The 
SOFC is modeled based on the following assumptions [28]: 

• Air is composed of 79% N 2 and 21% 0 2 , on a volume basis. 

• Temperatures at channel inlets are the same and, similarly, 
temperatures at the channel exits are the same. 

• Radiation heat transfer between gas channels and solid struc¬ 
tures is neglected. 

• The fuel cell is insulated perfectly, so there is no heat interaction 
with environment. 

• Only hydrogen reacts electrochemically, and CO is converted to 
C0 2 and H 2 via the water-gas shift reaction. 

• The gas mixture at the exit of the fuel channel reaches chemical 
equilibrium. 

• Contact resistances are negligible. 

• Unreacted gases are assumed to be fully oxidized in the 
afterburner. 

The GAX produces chilled water, and the generator is driven by 
exhaust gases from SOFC system. The energy and exergy analyses 
of the GAX are similar to those reported by Yari et al. [29], provid¬ 
ing a basis for validation of the present GAX analysis. The assump¬ 
tions used in the GAX model follow [29]: 
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Fig. 1. Proposed trigeneration system. 
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• The refrigerant is saturated at the channel exits of the 
condenser and the evaporator. 

• The solutions at the channel exits of the absorber and the gen¬ 
erator are at equilibrium and they are at the corresponding 
device temperatures. 

• Pressure drops due to friction in heat exchangers and in the 
connecting pipelines are neglected. 

• The approach temperature at either end of the GAX heat 
exchanger is assumed to be 0 K. 

• The system rejects heat to cooling water at the condenser and 
absorber. 


4. Modeling and analysis 

In this section, the modeling of the SOFC system, and the energy 
and exergy analyses of the considered trigeneration system, are 
described. The Engineering Equation Software (EES) is applied 
solve the resulting equations. 


4.1. Input data 

The input data applied in the analysis are listed in Table 1. 


4.2.1. Reformer model 

The ability of an SOFC system to simultaneously utilize hydro¬ 
gen and carbon monoxide as fuel is an important characteristic. 
Through direct internal reforming inside the fuel cell, carbon mon¬ 
oxide and methane result as a fuel mixture. The reactions that 
occur in the internal reformer are mainly endothermic, and attain 
the needed heat from the fuel cell. The use of an internal reformer 
reduces the dependence of the fuel cell on a cooling system. The 
following chemical reactions occur in the internal reformer [4]: 

CH 4 + H 2 0 — CO + 3H 2 (Reforming) (1) 

CO + H 2 0 -T C0 2 + H 2 (Shifting) (2) 

Eq. (1) depicts the steam reforming in which methane is converted 
to hydrogen. Eq. (2) presents the water-gas shift reactions. The 
hydrogen obtained from these reactions is used in the following 
relation electrochemical reaction the fuel cell: 

H 2 + ^0 2 —> H 2 0 (overall electrochemical reaction) (3) 

The mechanisms of the chemical and electrochemical reactions 
occurring at the anode and the cathode of the fuel cell are based 
on the following equilibrium equations [4]: 

x r -> [CH 4 + H 2 0 -> CO + 3H 2 ] (Reforming) (4) 


4.2. SOFC modeling 

For a methane fueled SOFC system, either an internal or 
external reformer is used. An internal reformer is usually superior 
because of its relatively low cost. The SOFC model is presented in 
the following two subsections. 


y r —> [CO + H 2 0 <-+ C0 2 + H 2 ] (shifting) 


( 5 ) 


Zr 


H2 + 2 ^2 


(Overall electrochemical reaction) (6) 


The molar conversion rates for Eqs. (4) and (5) are x r , y r , and z r , 
respectively. 


Table 1 

Input data. 


SOFC system [28] 


GAX system [29] 

Heating process heat exchanger (Heat recovery) 
Ambient condition 

Standard chemical exergy Szargut [30] 


Temperature difference between SOFC inlet and exit 

100K 

Fuel utilization factor 

0.85 

Active surface area 

0.01 m 2 

Baseline current density 

8000 A/m 2 

DC-AC inverter efficiency 

0.97 

Base inlet temperature to SOFC 

1000 K 

Exchange current density of anode 

6500 A/m 2 

Exchange current density of cathode 

2500 A/m 2 

Effective gaseous diffusivity through anode 

0.2 x 10 -4 m 2 / 

Effective gaseous diffusivity through cathode 

0.05 x 10" 4 m : 

Thickness of anode 

0.05 x 10 —2 m 

Thickness of cathode 

0.005 x 10 —2 n 

Thickness of electrolyte 

0.001 x 10 —2 n 

Thickness of interconnect 

0.3 x 10 —2 m 

Fuel compressor isentropic efficiency 

0.85 

Air compressor isentropic efficiency 

0.85 

Pump isentropic efficiency 

0.85 

Number of cells 

11,000 

Steam-to-carbon ratio 

2.5 

Afterburner combustion efficiency 

99% 

Fuel cell pressure drop 

2% 

Afterburner pressure drop 

3% 

Pressure drops for heat exchangers 

2% 

Efficiency of solution pump 

50% 

Effectiveness of heat exchanger 

80% 

Degassing value 

0.3 

Effectiveness of heat exchanger 

75% 

Pressure drop 

3% 

Ambient temperature 

298.15 K 

Ambient pressure 

101.3 kPa 

ch 4 

831.6 kj/mol 

co 2 

19.87 kj/mol 

CO 

275.1 kj/mol 

n 2 

0.72 kj/mol 

0 2 

3.97 kj/mol 

h 2 

236.1 kj/mol 

h 2 o (g) 

9.5 kj/mol 
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The fuel utilization factor and the air utilization ratio respec¬ 
tively can be defined as 


Up 

Uo 


( Fuel ■'conSLmm-d 

(Fuel) SU ppj e d 

(Air) consumed _ 

(Air)suppied 


_ i ^^iconsuined 
(^2) supplied 
(^(consumed _ 
(^2) supplied 


Zr 

3x r + y r 
Zr/2 
no 2 ,3 


(7) 

( 8 ) 


By applying mass balances to Eqs. (4)-(6), the molar flow rates n of 
the flowing gases can be determined as follows: 


n CH 4 ,ll — x r 


liman — r sc ■ iicH4,n — 2.5x r 


nH 2 ,i 2 = 3x r +y r -z r 


iico,i 2 = x r - y r 


tico 2 .i2 =y r 


tiH2o,i2 — ii H2o,n — x r — y r + z r — 1,5x r — y r + z r 


n 0 2 ,3 = 


Zr 

2Uo 



79 . 

Hn 2 , 3 = 21 x no 2 ,s 


tlN 2 .4 = tlN 2 ,3 


(9) 


tl3 = lto 2 ,3 + tlN 2 ,3 


1 14 = Ho 2 ,4 + Hn 2 .4 


n ll — n H 2 0,11 + ltCH 4 .ll 


tll 2 = ttH 2 ,12 + Hh 2 0,12 + tIC 0,12 + Hco 2 ,12 

The constants x r and y r are calculated using the equilibrium 
constant and current relations. The equilibrium constant for shift 
reaction can be written as 


ln/Cs = 


_ A g° 

R X Tp c, e 


In 


y,. x (3x r +y r -z r ) 
(x r -y r ) x (1.5x r -y r + z r ) 


( 10 ) 


where R, T FCi e and Ag° are the universal gas constant (8.314 
J/[mol K]), the temperature at the exit of the SOFC and the Gibbs 
free energy, respectively. The latter quantity for the shift reaction 
can be written as 


Aff° — pr° _]_ nr° _CT° _pr° 

^S s — Ss,C0 2 ^ &s,H 2 ®s,H 2 0 Ss,CO 

where 


( 11 ) 


g° = h-T F c,eS° ( 12 ) 

In Eq. (12), h and s° are the enthalpy per molar flow rate and the 
standard entropy per molar flow rate respectively. For ideal gases, 
h is function of temperature and s° is a function of temperature 
and standard pressure (101.3 kPa). 

The resulting current / and current density j respectively can be 
evaluated as follows: 


1 = J • A a 

. _ 2 ■ F ■ z r 


where A a , F and N FC denote the active surface area, the Faraday 
constant (96,485 C/[g mol]) and the number of cells. 

The work rate produced by the SOFC stack W ¥C , stack can be 
expressed as 

Wpctack = N fc ■ / ■ V c (15) 

where V c is cell voltage, defined as 

W = V N -V l0SS (16) 

Here, V N is the Nernst voltage and V^s the voltage loss, which is the 
sum of three separate voltage losses (activation, ohmic, and concen¬ 
tration). That is, the voltage loss can be written as 

V loss = V ohm + 1/act + V cone (17) 

The calculations of the Nernst voltage and the voltage losses are 
described in Table 2. There, p , L, R c , j oa , j oc , D ceff and D ae ff denote 
electrical resistivity of a cell component, thickness of a cell compo¬ 
nent, resistivity contact, exchange current density of anode, 
exchange current density of cathode, effective gaseous diffusivity 
through the cathode and effective gaseous diffusivity through the 
anode, respectively. Also, the subscripts a, c, e, int, cone, act and 
ohm denote anode, cathode, electrolyte, interconnect, concentra¬ 
tion, activation and ohmic, respectively. 

4.2.2. Validation of SOFC model 

The experimental data with methane fuel as introduced by Tao 
et al. [34] are used to validate the proposed model. The validation 
is depicted in Table 3, which shows the changes of both cell voltage 
and power density with the current density. The model is seen to 
exhibit good agreement with the experimental data, with an error 
of less than 6%. 

4.3. Energy analysis of system 

At steady state and neglecting changes in kinetic and potential 
energy, an energy rate balance for a system component can be 
expressed as 

Q-W my^hjhi - Yri e h e (18) 

i e 


Table 2 

Electrochemical equations. 


Type 

Equation 

Nernst potential 

^ = -«+% ? ln 

a H2 o = 

Ph 2 , 12 yH 2 ,12 ^12 Po 2 , 4 y 0 2 , 12 '^4 

«H 2 .12 - -4- - . “0..4 - 

A§ — Sh 2 o ~ Sh 2 — 2 %g 2 
g° = h-T FCs s° 

Ohmic loss [31] 

Vohm = (Rc + EiPih) i 

p„ = ( 3.34 X 10 4 exp(-10300/T FCe )^ 1 

Pa = (95 X 10 6 /Tfl;.e exp( —1150/Tfc.e)) 1 
p c = (42 X 10 6 /T FC , e exp(—1200/T Fc .e)j 1 

Pin, = ( 9.3 X io 6 /r fC , e exp(-noo/T FC , e )) _1 

Activation polarization [32] 

^7 act = V a ct,a + ^7 act,c 

V act . a =^- (sinh- 1 ^)) 


V^=5^-(sinh- 1 ^)) 

Concentration loss [33] 

V cone = V cone,a + ^7 cone,c 


V co „c,a=^.(ln(t +££)- In (l-i)) 


V conc.c = - (tP ' ln (! - 1) ) 


Jos = 2.F.Ph 2 .D ae ff/R.T FC . e ia 


Jcs = 4.F.P 02 .D ceJr /(( P l^t)R.T FC ,.L c ) 


(13) 

(14) 
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Table 3 

Comparison of present model results with results Tao et al. [34], 
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Current density (A/m 2 ) 

Cell voltage (V/) (Present work) 

Cell voltage (V0 (Tao et al.) 

Power density (W/m 2 ) (Present work) 

Power density (W/m 2 ) (Tao et al.) 

0.2 

0.79 

0.76 

0.158 

0.15 

0.3 

0.711 

0.68 

0.216 

0.21 

0.4 

0.644 

0.62 

0.253 

0.26 

0.5 

0.56 

0.57 

0.288 

0.295 

0.6 

0.51 

0.52 

0.3 

0.315 


where Q, W and h denote heat transfer rate, work transfer rate and 
molar enthalpy, respectively. This equation states that the total rate 
at which energy enters the control volume equals the total rate at 
which it exits. The expressions applied to determine the energy effi¬ 
ciency of the system are presented in Table 4. 


4.4. Exergy analysis of system 


Exergy is defined as the maximum work obtainable from the 
combined system of a system plus its environment, as the system 
passes from a given state to the dead state while interacting with 
the environment only. Alternatively, exergy can be viewed as the 
minimum work input required to bring the system from the dead 
state to a given state. Unlike energy, exergy is not conserved but is 
destroyed by irreversibilities. The exergy associated with a mate¬ 
rial flow is the sum of its physical and the chemical exergy. That is, 

Ex = Ex ph + Ex ch (19) 

where 


Ex ph = Jjii((hi - h 0 ) - T 0 (Si - s 0 )) 

i 

Ex ch = h (J>ex) h ’° + RT.^yJny. 


( 20 ) 

( 21 ) 


In Eq. (20), s is molar entropy and in Eq. (21), exf 10 is the standard 
chemical exergy of a species, as presented by Szargut [30] and listed 

in Table 1. 

The exergy rate balance of a control volume at steady state can 
be written as 


Table 4 

Energy efficiency definitions and system relations. 


Input energy 

Qin = 

h(TH4,ll ’ LHVcH4 

Stack AC power of fuel cell 

W FC,stack.ac — tfinv ' ^FC,stack 

Net electrical power of system 

w„ et 

— W FC,stack,ac ~ ( Wwp,/ + Wac 


+w fc 

+ Wmainpump + Wwp,//) 

Fuel cell efficiency 

*1 FC = 

WFCstack.ac-(W ac +Wfc+Wwpj) 1 QQ 
Qin 

Net electrical efficiency of system 

4el = 


Heating cogeneration efficiency 

^ lcog,h 

_ VVnet+Qb 1Q0 

Qin 

Cooling cogeneration efficiency 

tlcogt 

- . ioo 

Qin 

Trigeneration efficiency 

tftri — 

Wnet+Qeva+Qb Q0 

Electrical to heating ratio 

t el,h = 

= W„ et /Q h 

Electrical to cooling ratio 

te/,c = 

Wnef/Qeva 


Table 5 

Exergy efficiency definitions and relations for system. 

Input exergy 

EXj n - 

■ rr—c/i ,0 

1 HCH4,11 • ex CH 4 

Exergy efficiency of SOFC system 

'I'FC = 

(W FCsl ack,ac-(W wp j+W fc +Wac) m 

Net electrical exergy efficiency 

'Pel = 

if- 10 ° 

Heating cogeneration exergy efficiency 

^ cog,h 

_ lVnet+(EX47-E*46) -j QQ 

Cooling cogeneration exergy efficiency 

*A cog,c 

_ W ne t+(£X4i-£X4o) -J QQ 

Trigeneration exergy efficiency 

iptri = 

W net + (EX4 7 -Ex46 )+(Ex 4 , -EX 40 ) | QQ 


5.1. Current density 


i='E(l-^Qj-W cv + '£ E x i -'£Ex e ( 22 ) 

where i denotes inlet and e denotes exit. Also, Qj represents the rate 
of heat transfer at the location on the boundary where the instanta¬ 
neous temperature is T). The accompanying exergy transfer rate can 
be evaluated as (1 - T 0 /Tj)Qj. The term Ex, accounts for the exergy 
transfer rate accompanying mass flow and flow work at inlet i. Sim¬ 
ilarly, £x e accounts for the exergy transfer rate accompanying mass 
flow and flow work at exit e. The flow exergy rates Ex, and Ex e are 
evaluated using Eqs. (21) and (22). In Eq. (22), one-dimensional 
flow is assumed at locations where mass enters and exits. Finally, 
the term / accounts for the exergy destruction rate due to irreversi¬ 
bilities within the control volume. In Eqs. (20)-(22), the subscript 0 
denotes the value of a property at the surrounding conditions. The 
exergy efficiency relations for the system are listed in Table 5. 

5. Results and discussion 

The thermodynamic performance of the trigeneration system 
and the influence of variations of design parameters are examined. 
Since current density and inlet temperature to the SOFC signifi¬ 
cantly affect system performance parameters (e.g., energy and 
exergy efficiencies), we focus on them here. 


The current density is an important parameter in cycle 
performance. Figs. 2-8 show the influence of current density on 
the voltage, efficiency, power output, electrical to heating and cool¬ 
ing ratios, methane inlet flow rate, exergy efficiency and exergy 
destruction rate, for a constant SOFC inlet flow temperature of 
1000 K. In these figures, current densities from 5500 to 9500 
A/m 2 are examined. 

Fig. 2 indicates that the cell voltage V c decreases and the total 
voltage loss \/ loss increases as the current density increases. The 
total voltage loss is the sum of three voltage losses (activation 
polarization, ohmic and concentration), and the voltage losses are 
seen in Fig. 2 to be chiefly due to the activation polarization voltage 
loss. 

Fig. 3 shows that the efficiency of the SOFC decreases as current 
density increases, mainly because the SOFC cell voltage decreases 
with increasing current density. Also, the net power efficiency is 
observed in this figure to decrease as current density increases. 
The maximum net electrical cycle efficiency is 46% at 5500 A/m 2 
whereas the minimum is 31% at 9500 A/m 2 . The highest cooling 
cogeneration (integration of SOFC with GAX) efficiency is 57% 
while the lowest is 48%. Conversely, the maximum heating cogen¬ 
eration (integration of SOFC with heating load) efficiency, 67%, 
occurs at a current density of 5500 A/m 2 and the minimum, 60%, 
at a current density of 9500 A/m 2 . Fig. 3 reveals that the highest 
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Fig. 2. Variation with current density of several fuel cell voltages at T, = 1000 K. 

and lowest trigeneration efficiencies are approximately 79% and 
77%, respectively. It is also found that smallest gain in the effi¬ 
ciency of the trigeneration cycle compared with the net electrical 
efficiency is a considerable, at 33%. 

The influence of the current density is shown in Fig. 4 of the 
powers produced by the fuel cell and consumed by the pump as 
well as by the compressors. For a constant fuel utilization factor 
and SOFC inlet temperature, increasing current density is associ¬ 
ated with a linear increase of the fuel inlet molar flow rate (as 
shown in Fig. 5). That is, more chemical energy is converted to 
electrical energy and more cells current is produced. As the cell 
current increases and voltage decreases, a maximum SOFC-AC 
power output rate is thus obtained. It can be seen in Fig. 4 that, 
as current density increases, the net electrical power increases to 
a point and then decreases. There exists an optimal current density 
at which the maximum net electrical power is attained and, here, 
the maximum net electrical power is 421 kW at 8500 A/m 2 . The 
electrical power use of the air compressor rises with current den¬ 
sity due to increasing air flow rate through the compressor. 

Fig. 6 illustrates the influence of the current density on the elec¬ 
trical to heating and cooling ratios. The electrical to heating ratio is 
seen to range from almost 4 at 5500 A/m 2 to 1.1 at 9500 A/m 2 , 
while the electrical to cooling ratio varies from almost 4 at 
5500 A/m 2 to 1.7 at 9500 A/m 2 . 

The effect of current density on exergy efficiency is shown in 
Fig. 7. The net electrical exergy efficiency is observed to decrease 
with increasing current density, primarily because cell voltage 
decreases and therefore power output from the SOFC decreases 
as the current density increases, as observed in Figs. 1 and 3. The 
net electrical exergy efficiency is declines from 44% to 29%, for 
the range of current densities considered. The cooling cogeneration 
exergy efficiency is around 0.5% higher than the net electrical 
exergy efficiency, mainly due to the small flows of cooling energy 
compared to electrical energy in the system. Also, the heating 
cogeneration exergy efficiency is seen in Fig. 7 to decline with 
increasing current density, and the maximum and minimum heat¬ 
ing cogeneration exergy efficiencies are almost 46% and 32%, 
respectively. The trigeneration exergy efficiency is 0.5% or less 
higher than the heating cogeneration exergy efficiency. This effi¬ 
ciency increase is attributable to the small value of the cooling 
cogeneration exergy efficiency compared with the net electrical 
exergy efficiency. 

The variation with current density of the exergy destruction 
rates in the main system components are given in Fig. 8, where 
all component exergy destruction rates are seen to rise with 
current density. The exergy destruction rates of the air compressor, 
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Fig. 3. Variation with current density of several system energy efficiencies at 
Ti= 10001 <. 


the heat exchanger for the heating process, the absorber and 
desorber, the evaporator, the condenser and the fuel heat exchan¬ 
ger increase slightly (by less than 30 kW in total) as current density 
rises. But the variation in exergy destruction rate with current den¬ 
sity is significantly more notable for other components, with the 
highest exergy destruction rates in the air heat exchanger and 
the SOFC. The maximum exergy destruction rate of the SOFC is 
177 kW at 9500 A/m 2 and the minimum is 71 kW at 5500 A/m 2 , 
while the exergy destruction of the air heat exchanger changes 
almost from 114 kW at 5500 A/m 2 to 302 kW at 9500 A/m 2 . 

5.2. SOFC inlet temperature 

The effect of varying the SOFC inlet temperature on the perfor¬ 
mance of the trigeneration system is illustrated in Figs. 9-14, and 
is seen to differ from the effect of varying current density. The 
influence of SOFC inlet temperature is investigated here for a con¬ 
stant current density of 8000 A/m 2 . Fig. 9 shows that, as SOFC inlet 
temperature increases, activation and concentration voltages 
increase but the ohmic voltage decreases, so the cell voltage 
increases to a point and the decreases. There exists an optimal 
SOFC inlet temperature at which a maximum cell voltage is 
attained and here the highest cell voltage is almost 0.59 V at 
633 °C. 

Fig. 10 indicates that as SOFC inlet temperature increases, the 
energy efficiency of the SOFC increases to a point and then 
decreases. This trend mimics that of the cell voltage with varying 
SOFC inlet temperature (see Fig. 9). The highest net electrical cycle 
efficiency is 41.0% at 633 °C while the lowest is 26.8% at 500 °C. The 
maximum cooling cogeneration (integration of SOFC with GAX) 
efficiency is 55.0% and the minimum is 47.0. Correspondingly, 
the highest and lowest heating cogeneration (integration of SOFC 
with heating load) efficiencies, 64.3% and 56.5%, are obtained at 
SOFC inlet temperatures of 633 °C and 500 °C, respectively. The 
maximum trigeneration efficiency is seen in Fig. 10 to be 78.3%. 

The effect of varying SOFC inlet temperature on the powers pro¬ 
duced by the fuel cell and consumed by the pump as well as by the 
compressors is depicted in Fig. 11. With increasing inlet tempera¬ 
ture, the SOFC-AC output power increases to a maximum of 
503 kW at 633 °C and decreases above this temperature. The rea¬ 
son for this trend is similar to that discussed above regarding the 
influence of varying current density on cell voltage. The maximum 
net electrical power is almost 450 kW at 633 °C. 
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Fig. 5. Variation with current density of the methane inlet rate at T, = 1000 K. 
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Fig. 8. Variation with current density of system and component exergy destruction 
rates at T, = 1000 K. 
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Fig. 6. Variation with current density of electrical to heating and cooling ratios at 
T, = 1000 K. 


Fig. 9. Variation with SOFC inlet temperature of several voltages at j = 8000 A/m 2 . 





















326 


F. Ranjbar et al./Energy Conversion and Management 87 (2014) 318-327 



Inlet temperature to the SOFC (°C) 


Fig. 10. Variation with SOFC inlet temperature of system and component energy 
efficiencies at j = 8000 A/m 2 . 
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Fig. 11. Variation with SOFC inlet temperature of several power interactions at 
j = 8000 A/m 2 . 
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Fig. 12. Variation with SOFC inlet temperature of electrical to heating and cooling 
ratios at j = 8000 A/m 2 . 



Fig. 13. Variation with SOFC inlet temperature of system and component exergy 
efficiencies at j = 8000 A/m 2 . 



Fig. 14. Variation with SOFC inlet temperature of system and component exergy 
destruction rates at j = 8000 A/m 2 . 

The effects of varying SOFC inlet temperature on the electrical 
to heating and cooling ratios are shown in Fig. 12. These variations 
exhibit similar trends, and the maximum electrical to heating and 
cooling ratios are observed to be 1.76 and 2.92, respectively. 

The influence of SOFC inlet temperature on exergy efficiencies is 
shown in Fig. 13, where the trends exhibited for all cases shown 
are seen to be similar. Increasing SOFC inlet temperature t causes 
the exergy efficiencies to increase to a point and then decline. 
The maximum and minimum trigeneration exergy efficiencies 
are 42.6% and 28.7%, respectively. 

The effect of varying SOFC inlet temperature on the exergy 
destruction rates of the main system components are illustrated 
in Fig. 14, where the maximum exergy destruction rates are seen 
to occur in the air heat exchanger and the SOFC. As SOFC inlet tem¬ 
perature increases, the exergy destruction of the SOFC decreases to 
a minimum value of 126.7 kW at 667 °C and then increases 
slightly. 


6. Conclusions 

The energy and exergy assessments performed of a novel tri¬ 
generation system based on solid oxide fuel cell (SOFC) under 
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steady-state operation and using a zero-dimensional approach 
provide useful insights. The influence is characterized of several 
design parameters such as SOFC current density, SOFC inlet flow 
temperature and fuel utilization factor on numerous performance 
parameters. Several significant conclusions can be drawn from 
results of this research: 

• The energy analysis of the system demonstrates that trigenera¬ 
tion increases energy efficiency by at least 33% compared to a 
simple SOFC power cycle. That is, the maximum energy effi¬ 
ciency is 79% for the trigeneration system, 69% for heating 
cogeneration, 58% for cooling cogeneration and 46% for stand 
alone electricity generation. 

• The highest trigeneration exergy efficiency is almost 47% for the 
given conditions, and the exergy efficiencies decrease for the 
power cycle, cooling cogeneration, heating cogeneration and 
trigeneration as the SOFC current density increases. Further, 
the trigeneration energy and exergy efficiencies decrease with 
current density, and the maximum net electrical power is 
achieved at an optimal current density. As SOFC inlet flow 
temperature increases, the trigeneration energy and exergy effi¬ 
ciencies and net electrical power increase to a maximum and 
then decrease. 

• The main exergy destructions, which merit attention in efforts 
to improve system efficiency, occur in the air heat exchanger, 
the SOFC and the afterburner. Current density and SOFC inlet 
temperature significantly affect the exergy destruction rates in 
several of the main system components. 
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